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ABSTRACT 

The abundance of polycyclic aromatic hydrocarbons (PAHs) in low- and high- 
metallicity galaxies has been widely discussed since the time when detailed infrared 
data for extragalactic objects were first obtained. On the scales of entire galaxies, a 
smaller PAH abundance in lower-metallicity galaxies is often observed. We study this 
relationship for star-forming regions in nearby galaxies, for a sample containing more 
than 200 HII complexes, using spatially-resolved observations from the Herschel Space 
Observatory and Spitzer Space Telescope. We use a model for the dust emission to esti- 
mate the physical parameters (PAH abundance, metallicity, ultraviolet radiation field, 
etc.) of these complexes. The same correlation of PAH abundance with metallicity, as 
seen for entire galaxies, is apparently preserved at smaller scales, at least when the 
Kobulnicky & Kewley metallicity calibration is used. We discuss possible reasons for 
this correlation, noting that traces of less-effective PAH formation in low-metallicity 
AGB stars should be smeared out by radial mixing in galactic disks. Effective destruc- 
tion by the harder and more intensive ultraviolet field in low-metallicity environments 
is qualitatively consistent with our data, as the ultraviolet field intensity, derived from 
the infrared photometry, is indeed smaller in HH complexes with lower metallicity. 

Key words: infrared: galaxies ~ ISM, galaxies: photometry 



>. ' 1 INTRODUCTION 



Polycyclic aromatic hydrocarbons (PAHs) are now a well- 
recognized component of the interstellar medium (ISM), 
both in our Galaxy, and in external galaxies. Along with 
very small grains (VSGs, r ~ 100 A), PAHs are believed 
to be the dominant source of non-stellar emission in the 
near-infrared and mid-infrared ranges (2 — 20 /xm) in star- 
forming galaxies. This emission often constitu tes a sizable 
fract ion of a galaxy's total infrared luminosity (|Smith et al.l 
l2007h . and cannot be explained by large grains alone, which 
are too cold in the diffuse ISM to emit at these wavelengths. 
On the other hand, grains with radii much smaller than 
the typically- assumed value of 0.1 ^m can heat up to very 
high temperatures upon a single photon absorption. Thus, 
stochastically-heated VSGs can produce mid-infrared con- 
tinuum emission, while radiatively-excited bending and vi- 
brational modes in PAHs are the source of several spectral 
features in the ~ 2 — 20 /xm wavelength range. 

As an ultraviolet (UV) photon is needed to excite a 
PAH or VSG particle, the emission from t hese particles 
should be related to the overall UV emission l|Draine fc Li I 
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l200ll : lLi fc Draine |[20o3 ). Because of this expected connec- 
tion between the UV luminosity and the infrared (IR) lu- 
minosity, specific PAH bands, mid-IR continuum, and their 
combinations with other observables are now considered as 



possible indicators of star formation (e.g.. ICalzetti et al 
20081: iKennicutt et al.l 12003 : iTrever et al 



2007: Zhu et al 



20101 : iKim et al 



mm- 



However, the UV radiation that ex- 
cites PAH molecules can also destroy these particles, de- 
pending on the strength and hardness of the radiation field. 
The harsher UV fields expected in star-forming regions can 
therefore cause a significant local depletion of small parti- 
cles; thus, the relation between PAH emission and the star 
formation rate (SFR) can be more complicated than a sim- 
ple monotonic function. To assess the applicability of PAH 
emission as a tracer of the SFR, it is necessary to consider 
processes of both their formation and destruction in detail. 

Despite the significant observational progress in this 
field, some key theoretical questions about the PAH lifecy- 
cle are still debated. Of particular interest, both for PAH 
physics and for their usage as a SFR indicator, is the 
empirically-derived relation between th e metallicity of a 
galaxy an d the relative PAH abundance llEngelbracht et al 



2005; Madden et all 



20071 : IWu et al.l 12007 



20061: iDraine et all IIoOtI : ISmith et al 



Hunt et al.ll2010l ). Specifically there 
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appears to exist a metallicity threshold (12 + log(0/H) from 
~ 8.0 to 8.4 (depending on the adopted metaUicity calibra- 
tion) which separates galaxies into groups with high and low 
PAH abundance (in terms of the PAH fraction of the total 
dust mass, i.e. Qpah). 

Two general scenarios has been proposed to explain 
the origin of this threshold. If PAHs are synthesized in 
the atmosphe res of carbo n-rich asymptotic giant branch 
(AGB) stars l|Latterl [l99lf ). the efficiency of their produc- 
tion should depend on the carbon-to-oxygen ratio, which 
can be a function of metallicity. According to this scenario, 
PAHs are formed l ess efficiently by evo lved stars in low- 
metallicity galaxies (|Galliano et al.ll2008l ). According to the 
"destructive" scenario, PAH destruction can be more inten- 
sive i n low-metallicit y environments due to the harder UV 
field J^-"^ -1'"^ 



Galliano et al.'''2005'; 'Mad den et"aLll2006l ) or stronger 



shocks ([O'Halloran et al. 2006). As neither of these scenarios 
now seems clearly preferable, the question of the formation 
and destruction of large molecules in these environments re- 
quires additional study. 

Studies of PAHs have been tightly related to the devel- 
opment of infrared telescopes. Our knowledge about their 
presence and importance has expanded significantly, thanks 
to the Infrared Space Observatory (ISO) and especially to 
the Spitzer Space Telescope. The Spitzer IRAC and MIPS 
passbands cover PAH spectral features, in addition to con- 
tinuum emission from larger grains. The combination of ob- 
servations at short (3.6 — 24 fim) and long (70 and 160 /xm) 
wavelengths makes it possible to determine the PAH abun- 
dance relative_tothB_total dust mass in extragalactic ob- 
jects (jPraine fc Lill2007l ). However, for many galaxies, only 
unresolved photometry can be performed with Spitzer at 
longer wavelengths, where the spatial resolution (~ 38") 
is not sufficient to measure individual fluxes from smaller 
regions. It is thus not possible to resolve individual HII re- 
gions and star-forming complexes at these wavelengths in 
most galaxie s, except the for Large and Small Magel lanic 
Clouds (e.g.. ISlater et al. 2011 : Sandstrom et al.ll20ld . cor- 
respo ndingly), MlOl jOprdon et al.ll2008l ). and some others 



fe.g. .lHavnes et al.ll201ol) 



iMufioz-Mateos et al.l (|2009t ) performed a spatially re- 
solved s tudy of the dust prope rties of galaxies in the SINGS 
sample l|Kennicutt et al.]|2003l ) on a scale larger than that of 
individual HII complexes, using azimuthally-averaged emis- 
sion proflles. They found that the correlation between the 
metallicity and PAH content is observed not only in a galaxy 
as a whole, but also radially. They conclude that at large 
spatial scales this correlation is more readily explained by 
evolutionary effects. In that work, no distinction is made be- 
tween PAHs in star-forming regions and PAHs in the general 
diffuse ISM. 

We can expect that the details of PAH physics will 
be revealed with more clarity in spatially-resolved obser- 
vations of individual star- forming regions and complexes, 
rich in bright UV-sources, both heating PAHs and destroy- 
ing them. The intensity of 8 fim emission, and its correla- 
tion with metallicity, may depend on the source of heating. 
Within (or close to) HII regions, PAHs are heated by single 
ultraviolet photons produced by young stars. However, at 
some distance from these regions, the radiation field from 
older stellar populations can heat PAHs as w ell, and their 
emission is thus no longer a direct SFR tracer (jBoselli et al.l 



|2004 iBendo et al.ll20o4 ICrocker et al.ll2013l ). Therefore, it 
is desirable to separately observe PAH emission both in the 
vicinity of HII complexes, and in the general diffuse inter- 
stellar medium. 

To achieve this goal, we need the ability to study indi- 
vidual star-forming regions and complexes, both at the mid- 
and far-IR wavelengths. The resolution at 70 and 160 fim 
provided by the Herschel Space Observatory (about 5.2" and 
12", respectively) partly satisfies the requirements. Corre- 
sponding physical scale is about several hundred parsecs for 
nearby galaxies. This resolution in most cases does not allow 
studying individual HII regions but is sufficient to analyse 
large HII complexes and filaments. Thus, the Herschel data 
successfully complement the Spitzer data at shorter wave- 
lengths, which have a resolution of less than 1-2" for IRAC 
the passbands (3.6-8. 0/xm), and ~ 7" for MIPS 24 ^m. 

Previously, the spatially resolved analysis of the mid- 
IR proper ties and their r elatio n to the metallicity was per- 
formed bv lCordon et al.l (|2008l ) for a large sample of HII re- 
gions in MlOl. The goal of this paper is to extend this sample 
significantly and to study the properties of PAH emission in 
individual HII complexes within several nearby galaxies at 
various metallicities, using data from the Spitzer and Her- 
schel telescopes. We consider possible correlations between 
the various derived properties of the HII complexes, and dis- 
cuss the main scenarios for PAH formation and destruction 
in the context of these observations. 



2 OBSERVATIONAL DATA 

For this work, we consider 24 galaxies which were ob- 
served with both the Spitzer Space Telescope and the Her- 
schel Space Observatory. All of the Spitzer data exam- 
ined in this paper were obse rved as part of the SINGS 
survey (|Kennicutt et al.ll2003D . and most of the Herschel 
data were observed as part of the KINGFISH survey 
(|Kennicutt et al.]|201l[ ). although we also include data for 

3 galaxies (NGC 2403, NGC 6822, NGC 5194) which were 
observed with Herschel as parts of other programs. In Ta- 
ble[l] we list the name and morphological type of each galaxy 
(from Moustakas e t al.. .2010. ) , together with the proposal 
identification number (Prop. ID) and principal investigator 
for the Herschel data presented here. We also list the num- 
ber of HII complexes, analysed in each galaxy, and references 
to the works used to identify them. In the last column we 
indicate the physical resolution for each galaxy probed by 
the Herschel instruments. For most galaxies, the resolution 
is about a few hundred parsec, however, in some cases it is 
about 200 pc or even less which could mean that we probe 
individual giant HII regions. 

For the galaxies in Table [TJ we consider only data ob- 
tained with the IRAC and PACS (|Poglitsch et al.ll2010l) in- 
struments on Spitzer and Herschel, respectively. We down- 
loaded the reduced IRAC images at wavelengths of 3.6, 
4.5, 5.8 and 8.0 nm from the SINGS project websiteQ- For 
the Herschel data, we consider images at 70 and 160 fim, 
obtained with the PACS photometer. We downloaded the 
so-called Level 1 data from the Herschel Science Archival, 



^ http://sings.stsci.edu 
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which were reduced with version 6.1.1 of the automated re- 
duction pipeline. We produced final maps from the Level 
1 da ta using version 16 of the Scanamorphos algorithm 
l|Roussel.2012 ). with a pixel scale of 2'.'0. 

3 METHODS OF MEASUREMENTS AND ANALYSIS 

3.1 Photometry 

Since the point spread functions (PSFs) of the observations 
from different instruments and passbands differ from each 
other, we convolved all shorter-wavelength images to the 
resolution of 160 fim FAGS images, us ing the IDL convolu - 
tion procedure and kernels provided bv lAniano et al.ll|201lD . 
Finally, all the images were resampled to have a pixel scale 
of 2'.'0. 

The convolved and resampled images were used to per- 
form aperture photometry in the following fashion: the flux 
in an aperture is measured for each HII complex f or which 
the m etallicity is available in the work by Mousta kas et al.l 
The aperture radius is chosen separately for each HII 
complex, depending on its size, but they are not less than 
the size of the lowest resolution element (12"). Therefore, no 
aperture corrections were applied. The centre of an aperture 
is chosen geometrically. Finally, we correct the flux for the 
contribution of partial pixels. 

To subtract the background emission, we estimate its 
average value and standard deviation, (Tring, within a ring 
having a width of a few arcseconds around the aperture. To 
exclude bright sources (e.g., neighbouring HII complexes) 
from the ring, we take into account only pixels with inten- 
sities smaller than 3(Tring. As the shift of an aperture by 1 
pixel may significantly (by ~ 10%) change the estimated 
background, we move the centre of each aperture in all di- 
rections by 1 pixel (in steps of 0.5 pixel) and calculate the 
average flux for all aperture positions. This reduces the influ- 
ence of the position error, too. The standard deviation from 
this average is o-pos. The two uncertanties, cjpoa and (Tring, 
are comparable in magnitude and in most cases contribute 
more or less equally to the total uncertainty. The final error 
of the aperture flux, cr, is then equal to 

= \/'^?ing + '^poB- (1) 

This value varies from 1% to 10% for most complexes, but 
can be larger for faint complexes that cannot be clearly sepa- 
rated from neighbouring complexes. Locations with a signal- 
to-noise ratio < 3 at 8 fim are excluded from the study, which 
means that only large and bright complexes are analysed, as 
small and faint complexes are discarded because of the con- 
volution process. The results of photometry for all bands are 
presented in Table [2] (see below). 

Finally, we note that measuring reliable fluxes for low- 
metallicity galaxies, like DDO 053 and Holmberg I, is com- 
plicated by the extremely low signal-to-noise ratio. As a re- 
sult, we had to abandon most complexes from these galaxies, 
but we present fluxes for some of the brighter complexes in 
Table [2] (see below). 

3.2 Fitting HII Complex Spectra 

We apply the algorithm that has been used in lDraine et all 



(|2007l ) to calculate global SEDs of SINGS galaxies, but use it 
to compute emission spectra of individual HII c omplexes. To 
model a dust emission spectrum, we adopt the iDraine fc Lil 
( 200 7) approach but use an alternative method of comput- 
ing temperature distribution functions. Details are provided 
below. 

We assume that a theoretical emission spectrum of 
an HII complex, l ocated at the distance D, is given by 
(|Draine et al.ll2007l ) 

model = Q.*Bi,{Ti) -I- -^^^ [(1 ~ 7)PJ!('7PAH, Ui-niu) 

+ 7pi(9PAH, C/min)] • (2) 

Here the first term with T* = 5000 K describes the stellar 
contribution, Mdust is the total dust mass in the HII complex 
and S7t is the starlight dilution parameter. 

The radiation field in a complex is described by the 
scaling factor U , which is given in units of the interstel- 
lar radiation field in the solar vicinity, as estimated by 
iMathis et al. Thus, the energy density of the radi- 

ation field is U X For U, we use the representation 
given by Eq. 23 from iDraine fc Lil l|2007h . In this represen- 
tation, a fraction 1 — 7 of all the dust is exposed to the 
"minimum" radiation field C/min, while the remaining dust is 
illuminated by the enhanced radiation field, having a power- 
law distributio n with an exponent of a, and an upper limit 
of [/max. As in lDraine et al.l (|2007l ). we adopt the fixed val- 
ues [/max = 10® and a = 2. Functions pi and pi represent 
emission per unit frequency per unit mass of dust (with a 
PAH fraction qpah) illuminated by the minimum radiation 
field and by the enhanced radiati on field, resp e ctivel y. 

We adopt a grain model from IDraine fc Lill|2007l l which 
represents a dust mixture consisting of carbonaceous par- 
ticles (including VSGs and PAHs) and amorphous sili- 
cates. The Milky Way dust with _Rv = 3.1 is utilized. For 
grains larger than 250 A an eqilibrium temperature is cal- 
culated for a given starlight intensity. For smaller grains 
we take the stochastic heating into account and calcu- 
late the temperature distrib ution functions as described by 
iPavlvuchenkov et all l|2012l) . 

We fol l ow t he fitting procedure described in 
IDraine et all 1120071 ). The best- fit SED for an HII complex 
is found using the observed fiuxes at 3.6, 4.5, 5.8, 8.0, 
24, 70, and 160 /im. The fitting parameters are the PAH 
mass fraction qpah, the minimum starlight intensity [/min, 
the parameter 7 describing the fraction of dust headed by 
radiation field with intensity greater than f/min, the dust 
mass Mduat and the solid angle fi* of stellar emission. The 
last two parameters are used for normalization only. 

In the parameter grid, used for fitting, the gpAH val- 
ues range from 0.47 to 4.6% (a step is 0.1); the minimum 
starlight [/min varies from 0.1 to 25 (a step is 0.1); and 7 
varies from to 0.3 (we sample this parameter 10 times 
within this range in logarithmic scale). The best-fit model is 
the model with the minimum value of , where is given 
by 

/ rohs rmodcl \2 

2 X ^ l../i^,band 7i^,band/ /q\ 

^ ^ 2-^ ('„obs \2 I ('_modcl-i2 ' ^ > 

band band,* I" band ) 

Here /"Iband and (jjj!jJd are the observed fiux density at 
a given passband and the uncertainty of this value, while 
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Table 1. Sample of galaxies 

Object Prop. ID Proposer Type Number of Ref. Physical 

complexes resolution, kpc 





KPOT rkennini 1 


R. Kennicutt 


Im 


2 


5 





21 




KPOT rkpnnini 1 


R. Kennicutt 


lABm 


I 


5 





22 


HolmbGrg II 


KPOT_rkennicu_l 


R. Kennicutt 


Im 


Q 


20 





20 


IC 2574 


KPOT rVpnnirn 1 


R. Kennicutt 


SABm 


7 


5 21 





23 


NGC 628 


K Pf vlr^iTiTi 1 fii 1 


R. Kennicutt 


SAc 


ll 


1,2,3,4 





42 




KPOT rkpnnini 1 


R. Kennicutt 


SABd 


23 


2 


Q 


53 


NGC inqy 


S(T )P T" Vr^vi n ir'i 1 "-t 


R. Kennicutt 


SBb 


g 


5 





99 


Non 24ns 


KPCT rwilsofll 1 


C. Wilson 




ll 


7 8 





18 


NGC 3184 


KPOT_rkennicu_l 


R. 


Kennicutt 


SABcd 


14 


2,9 





65 


NGC 3198 


KPOT_rkennicu_l 


R. 


Kennicutt 


SBc 


8 


9 





80 


NGC 3351 


KPOT_rkennicu_l 


R. 


Kennicutt 


SBb 


6 


1,15 





54 


NGC 3521 


KPOT_rkennicu_l 


R. 


Kennicutt 


SABbc 


8 


1,9 





58 


NGC 3621 


KPOT_rkennicu_l 


R. 


Kennicutt 


SAd 


21 


9,16 





38 


NGC 4254 


KPOT_rkennicu_l 


R. 


Kennicutt 


SAc 


10 


3,17,18 





96 


NGC 4321 


KPOT_rkennicu_l 


R. 


Kennicutt 


SABbc 


7 


3,18 





83 


NGC 4559 


SDP_rkennicu_3 


R. 


Kennicutt 


SABcd 


12 


9 





60 


NGC 4725 


KPOT_rkennicu_l 


R. 


Kennicutt 


SABab pec 


7 


9 





69 


NGC 4736 


KPOT_rkennicu_l 


R. 


Kennicutt 


SAab 


7 


1,14 





30 


NGC 5055 


KPOT_rkennicu_l 


R. 


Kennicutt 


SAbc 


5 


3 





45 


NGC 5194 


KPGT_cwilso01.1 


C. Wilson 


SABbc pec 


14 


10,11 





44 


NGC 6822 


SDP_smaddc01.3 


S 


Madden 


IBm 


5 


12,13 





02 


NGC 6946 


KPOT_rkennicu_l 


R. Kennicutt 


SABcd 


7 


3,4 





40 


NGC 7331 


KPOT_rkennicu_l 


R. Kennicutt 


SAb 


2 


1,14 





84 


NGC 7793 


KPOT_rkennicu_l 


R. Kennicutt 


SAd 


9 


3,19 





23 



Refer ences: (11 iBresolin et a.1.1 lll999h . (2) Ivan Zee et al.l lll998h . ( 3) iMcCall et alJ lll985l) , (4 ) iFerguson et al l lll998h . (5 ) ICroxall et al. 
1 20091') . (6) iLee fc Sk illman' l'2004'l. (71 Garnett et al.' ('19991 '). (81 iGarnett et alF l(l997''). (9) Zaritsky et al.' ('1994'), (10) Brcsol in et al, 
( 2004 1. (i nlDiaz et al. (1991), (12) Lee et al (200 6), (1 31 Hodee et al.l l ll988| l. (14) Ocy & Kennicutt (1993), (15) Brcs olin & K ennicutt 
I2OO2I1. 16lRvderl (Il995l. (171 iHenrv et al.l |l994l l. (181 IShields I (|l99ll l. (191 lEdmunds fc Pagell (|l984j l. (201 iHodge et al.l (jl99i l. (211 
lArsenault fc Rovl (Il986i 1 



/™band '^bancT' the theoretical flux de nsity and its 
"unce rtainty," which has been proposed by IPraine et al.l 
(|2007h to account for different levels of uncertainties in dif- 
ferent bands. It is estimated as 0.1/™band- The best-fit 
values vary from 0.5 to 3, depending on the region. Uncer- 
tainties in the derived parameters are calculated using the 
Monte-Carlo method. 



3.3 Metallicities of the HIT complexes 

Many different methods have been developed to determine 
the metallicity of HII regions. One of the main methods is 
the so-called "direct" To method, based on the electron tem- 
perature measurement (e.g. Stasihska 2007). However, this 
method requires observations of very faint nebular lines, and 
consequently is of limited use. Besides this, various other 
"indirect" methods can be used. "While numerous methods 
are available, metallicity estimates obtained with different 
methods may differ significantly (up to 0.7 dex). Thus, re- 
gardless of which particular method is chosen, it is impor- 
tant for consistency to use the same method to evaluate the 
metallicities for all HII regions. 

We use the metalli cities of HII comp l exes in the SINGS 
galaxies calculated by iMoustakas et al] (|201(]| 1. They col- 
lected spectroscopic data from various sources and used two 
methods to estimate metallicities. These methods are the 
theore tical relationship presented by iKobulnickv fc KewlevI 
(|2004l l (hereafter, KK04) and the empirically-derived rela- 



tionship of jPilvugin fc ThuanI (120051 1 (hereafter, PT05). In 
both methods, relatively strong nebular lines ([Oii]A 3727; 
[Oiii]AA 4959, 5007; H^) are utihzed. 

These methods are convenient for a number of rea- 
sons. First, estimates obtained with them tend to bracket 
metallicities determined by various other methods. In other 
words, the "true" value of metallicity presumably lies some- 
where within the range bounded by values calculated with 
the KK04 and PT05 methods. Second, while other indirect 
methods are recommended for objects with specific values 
of luminosity, metallicity, etc., the KK04 and PT05 methods 
are more or less universal. The HII complexes in our sample 
possess a wide range of physical properties, so this "uni- 
versality" is desirable. For further discussions and a com- 
parison of these metho ds, we refer to the publication by 
IMoustakas et al.l (|2010l '). but some points are worth to be 
repeated here. 

First, metallicities determined with the KK04 method 
are systematically higher than those determined with the 
PT05 method (Fig. [TJ . As the latter metallici ties are similar 
to m etallicities obtained with the To method (|Egorov et al. I 
l2012h . the KK04 method results are also systematically 
higher than "direct" metallicity estimates. Second, when 
the PT05 metallicities are plotted versus galactocentric ra- 
dius, they show shallower grad ients and much more scat- 
ter than the KK04 metallicities (|Moustakas et al.|[201oh . We 
see the scatter in our result s, too (see below). Furthermore, 
iLopez-Sanchez et al.l (|2012l ) compared these methods, and 
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Figure 2. Comparison of gpAH (a) and !7niin (b) values obtained 
in this work and in Drainc et al. (2007). 



Figure 1. Comparison of metallicities estimated by KK04 and 
PT05 methods. A straight hne y = x is added for clarity. 



grain size distributions from IWeingartner" fc Draind (|200lD . 
and they provide parameters for distributions only for dust 
mixtures in this range of qpah- Thus, we are unable to mea- 
sure gpAH values lower than 0.5% and greater than 4.6%. 
Note that qpAH is insensitive to the band-to-band variations 
which can be caused by changes in the mean ionization state 
or grain size distribution. 



suggested that the KK04 method produces the most reli- 
able metallicities for stellar populations. 



4 RESULTS 

The results of the aperture photometry are presented in Ta- 
bled with the following columns: (1) sequential number; (2) 
designation of the HII complex; (3)- (4) right ascension and 
declination of the complex; (5) an aperture size; and (6)- 
(12) the fluxes at IRAC 3.6, 4.5, 5.8, 8.0, MIPS 24, PACS 
70, and PACS 160 \im bands. The results of fitting proce- 
dure are presented in Table |3l which contains the follow- 
ing columns: (1) sequential number; (2) designation; (3) the 
relative PAH mass abundance, Qpah; (4) 7 parameter; (5) 
minimum starlight intensity i 7min; (6) and (7) metall icities 
of the complex, determined bv lMoustakas et al.l (|201Cll ) with 
the KK04 and PT05 methods. 

The values of Qpah, 7 and C/min h ave been determine d 
for each galaxy in the SINGS sample bv lDraine erall l|2007f ). 
In Fig. [21 we compare the values of qpah and Umin for each 
galaxy, obtained by simple averaging of the values from 
Table [3l with the res pective values given by iDraine et all 
(|2007l ). Note that in I Praine et all (120071 ) the parameters 
were obtained based on global photometry of galaxies in- 
cluding diffuse medium not associated with HII complexes 
which can contribute a significant fraction to the total in- 
frared emission. The PAH abundances in both studies agree 
well with each other, while our starlig ht intensities are sys - 
tematically higher than the values of iDraine et all (|2007l ) 
(the same is true for 7 as well). This is expected, as we only 
consider HII complexes, and our Umin starlight intensities 
are biased toward higher values. 

Since we are mostly interested in PAHs, below we con- 
sider several relations between the metallicity and parame- 
ters related to PAHs. Note that complexes with minimum 
and maximum gpAH values represent, correspondingly, up- 
per and lower limits. In this study we use the models for 



4.1 Derived parameters of HII complexes 

We have determined mass fraction of PAHs in more than 
200 HII complexes for which metallicities are known. The 
combined results are presented in Fig. [3] both for the KK04 
metallicity calibration and for the PT05 metallicity calibra- 
tion. The metallicit y thr eshold inferred in previous stu dies 
(lEngelbracht et all boosi : fPraine et all l2007l : iGordon etHI 
|2008| ). which divides HII complexes with high and low PAH 
abundances, is apparently present in Fig. [3| It is ~ 8.4 in 
the KK04 scale and ~ 8.0 in the PT05 scale which is consis- 
tent with previous estimates. Depending on the method of 
metallicity estimation some pre vious studies quoted differ- 
ent values of the threshold. E. g., IDraine et~all (|2007l ). usmg 
the PT05 method , found the threshold to be ~ 8.1, while 
lEngelbracht et aL I 12005!) give the threshold value of about 
8.2, using empirical methods that are close to the PT05 one. 

When the KK04 calibration is used, a tendency is seen 
for the abundance of PAHs to grow with metallicity above 
the threshold (Fig. (3^). For the PT05 calibration there is 
more scatter on the gpAH-metallicity plot that has a stag- 
gered appearance (Fig. (SIj). As this scatter would blur any 
correlations that may exist in individual HII complexes, for 
the remainder of the present work, we adopt the KK04 
metallicities. It must be kept in mind, however, that our 
conclusions are based on the specific metallicity calibration. 

In Fig. (4] we show the metallicity vs. the PAH mass 
fraction qpah for several individual galaxies. In each subplot, 
we show the values for the entire sample (over all galaxies) 
as grey points, and the values for each individual galaxy 
(indicated by the caption) as black triangles. 

We have only a few low-metallicity HII complexes in our 
sample, partly because these complexes are too faint in the 
infrared range. Furthermore, the SINGS/KINGFISH sample 
includes only four low-metallicity galaxies, and none of the 
HII complexes in these galaxies cover the "critical" metallic- 
ity range from about 8.2 to 8.5 (in the KK04 scale), which 
apparently includes the turn-off point. However, at higher 
metallicities, the tendency for the PAH content to grow with 
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12 + log(0/H) is preserved, even when we consider individ- 
ual HII complexes instead of an entire galaxy. This can be 
seen in Fig. [4l where we show data for selected galaxies from 
our sample with metallicities above the "critical" range. De- 
spite the small number of data points for each galaxy, we still 
believe that our data are sufficient to conclude that the de- 
pendence of PAH content on the metallicity in star-forming 
galaxies is observed both globally and locally, at the scale of 
individual HII complexes. T his result is co nsistent with the 
similar conclusion made by [Gordon et al.l (,2008 ') for MlOl 
galaxy. 

In Fig. [S^ we show how the PAH mass fraction gpAH is 
related to the UV strength C/min. Solid circles correspond to 
HII complexes with gpAH less than 4.6%, which is the highest 
value we can measure with our method. Open circles are HII 
complexes where the plotted gpAH values represent lower 
limits, complexes with higher t/min tend to have lower gpAH, 
albeit with significant scatter. The only noticeable exception 
is the CCM91 region in M51, which has a lower limit for 
gpAH of 4.6%, and i7min = 21. It is interesting that high 
f^min values do indeed show some preference toward lower 
metallicity HII complexes (Fig. [Sja). The apparent grouping 
of points in the lower part of Figure simply reflects the 
lack of observed complexes and galaxies in the "critical" 
metallicity range. Probably, these intermediate systems are 
too faint in the infrared to become a part of some large-scale 
survey, but at the same time too metal-rich to deserve an 
individual study. 



4.2 Dust indices, PAH content, and metallicity 

Apart from the SED modelling, some indices can b e used 
to characterize the dust content. iDraine fc Lil (|2007l ) intro- 
duced the Pg.o and P24 indices 



Ps.o = 



iyF^{70fim) + !/F„(160^im) ' 



(4) 



(5) 



where -/^"^(S.O/xm) and F°''(24/im) are non-stellar fluxes at 
8 and 24 fim, correspondingly. Stellar emission signiflcantly 
affects Spitzer 3.6 /xm band; at 8 fim and 24 /xm, the stellar 
contribution is less important, but still should be taken int o 
account. We use the relations, given bv lHelou et al] l|2004h . 
to subtract stellar emission. At 70 and 160 [xm, the stellar 
contribution is negligible, and no correction is applied to 
( 70/xm) and F^(160/xm ) fluxe s . 
iMunoz-Mateos et al.l (|2009l ) proposed an analytical re- 
lation (their Eq. A4) between qfah, Ps.o and P24, which 
allows estimation of the PA H content in cases when gpAH 
goes beyond the limits of the lDraine fc Lil (|2007f ) spectra. In 
our data, the dependence of (/pah on infrared indices look 
similar, but our gpAH values are systematically lower than 
those computed for t he same Ps.o and P24 values with the 
analytical relation of iMunoz-Mateos et al.l (|2009l ). This is, 
probably, again related to the fact that we only consider HII 
complexes with enhanced UV radiation. In these complexes 
fewer PAHs are needed to produce a certain level of 8 //m 
emission compared to cases when this emission is averaged 
over the galaxy as a whole (or over a significant portion of 
the galaxy). 

We considered several relations between metallicity and 
the Pg.o and P24 indices. In Fig. (6] we show how Pg.o and 
P24 are related to the metallicity of an HII complex. Pg.o 
shows nearly the same correlation with metallicity as qpah 
(Fig. [6^). On the other hand, P24 appears to be totally un- 
correlated with metallicity (Fig. [SJa). 

Naturally, Pg.o is correlated with gpAH as well (Fig.[7K), 
but this index, unlike (?pah, is not bound to the model grid, 
so we used it to mark some gpAH values as lower and up- 
per limits. Specifically, we estimated the bordering values 
Pg"!]" and PgX that correspond to qpah of 0.5% and 4.6%, 
respectively. If Pg.o for an HII region is greater than Pg"o"' 
or less than Pg^o" we assume that the corresponding gpAH 
value is an upper or lower limit. 

The absence of correlation between P24 and metallic- 
ity implies that qualitative assessment of the PAH content 
in HII complexes can be done, using the Pg.o/P24 ratio or, 
equivalently, the Pg^oZ-^M ratio. This ratio is sometimes con- 
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Figure 4. Relation between gpAH and metallicity for individual HII complexes in selected galaxies. The grey points in each subplot 
show the values from the entire sample, while the triangles show only those values corresponding the galaxy indicated in the upper left. 
Galaxies which are not included in this plot have high metallicities so that the corresponding points are concentrated in the upper right 
part of the plot. 



sidered to be a quantitat ive measure of the PAH abundance 
l|Engelbracht et all lioosh . In our results, it is indeed well 
correlated with qpAu (Fig. [TJs). This greatly expands the 
ability to quantify the PAH evolution in HH complexes, as 
the Fg^Q I F2I ratio does not require photometry at longer IR 
wavelengths. 



5 DISCUSSION 

Our results indicate that the correlation between the PAH 
abundance and metallicity may exist both globally (for en- 
tire galaxies) and locally (for individual star-forming com- 
plexes). Here, we examine how the local character of the 
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Figure 6. Relation between Pg.o and P24 indices and the metallicity of an HII complex. 



qpAK — (0/H) correlation is related to our ability to under- 
stand the processes of PAH formation. 

According to the "AGB" scenario, a low gpAH value in 
low- metallicity environments reflects certain details of their 
synthesis in metal-poor AGB stars. Specifically, AGB stars 
populate low- metal licity galaxies as well as high-metallicity 
galaxies (|Tosill2009l ). but the carbon-to-oxygen ratio in the 
former may be too low for effective PAH formation. 

Neither observational nor theoretical data currently 
seem to support this assertion unequivocally. There are indi- 
cations that C/O ratio in the interstellar medium is roughly 
constant at metallicities lower than 12 -I- log(0/H) ~ 8.0 (at 
the direct m ethod scale) and g rows with (O/H) at higher 
metallicities (|Henrv et al.ll2000f ). This behaviour is reminis- 
cent of the observed PAH behaviour, and it is tempting to 
assume that the correlation between gpAH and (O/H) simply 
reflects the correlation between C/O and (O/H). 

However, a signiflcant number of extremely metal-poor 



stars with carbon overabundance implies that dependence of 
the C/O ratio on metallicity is more complicated than just 
"constant, then growing." Also, there is still no direct obser- 
vational evidence that PAHs are formed in the atmospheres 
of AGB stars, even though PAHs are observed in vicinities 
of some post-AGB and AGB stars with hot companions. 

Circumstantial evidence is more difficult to interpret, 
as there are at least two additional factors that make the 
situation confusing. First, observations indicate that acety- 
lene molecules, believed to be precursors for PAHs, are more 
abundant in low-metallic ity AGB stars llWoods et al. II2012I . 
and references therein). IWoods et al. I I 2012h have shown 
with chemical modelling that these molecules, and other 
hydrocarbons, form more efficiently in the atmospheres of 
metal-poor stars compared to stars of solar metallicity. Sec- 
ond, metal-poor stars not only produce abundant complex 
hydrocarbons, but ma y also expel them int o circumstellar 
space more efficiently (|Mattsson et al.ll2008l ). Also, from a 
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more general point of view, the evolution of C and O in 
AGB stars is defined by many processes that involve both 
synthesis and destruction, so that surface C/O ratio (sup- 
posedly related to PAH synthesis) depends not only on the 
star initial metallicity and mass, but also on tim e within 
the A GB phase (e.g., iLattanzio fc Forestini|[l999l : iHerwigl 
So, it is well expected that initial metallicity would 
not play a simple role in determining the PAH synthesis. At 
the very least, low metallicity by itself does not guarantee 
fewer PAHs. 

Another complication is the non-uniform and (possi- 
bly) time-dependent metallicity distribution in most late- 
type galaxies. If PAHs are synthesized in AGB stars with an 
efficiency that is an increasing function of metallicity, then 
the PAH content should depend on the metallicity of these 
parent stars (i.e., on the metallicity at locations where these 
stars were formed). The original relation between the PAH 
content and metallicity would be preserved if PAHs and/or 
their parent stars are not able to travel in the galaxy over 
significant distances (i.e., distances, over which the metal- 
licity changes appreciably). 

The question is if AGB stars can migrate over kiloparsec 
scales in their parent galaxies. Calculations show that radial 
mixing of stars can fiatten the metallicity gradient both in 
the Milky Way-type galaxies, and in less-massive galaxies, 
over a time span comparable to the lifeti me (from about 100 
Myr t o several Gyr) of an AGB star (ISellwood fc BinnevI 
l2002l : iMinchev et all I2OIII '). The metallicity range in the 
galaxies from our sample does not exceed a few tenths. To 
destroy the original correlation between the PAH content 
and the ISM metallicity, the star would have to have moved 
radially by a few kiloparsecs prior to the AGB phase, which 
seems probable according to theoretical models. If this were 
the case, we would expect the metallicity gradient of evolved 
stars to be shallower than the gradient of the ISM metallic- 
ity, or the gradient measured with young objects like O stars 
or open clusters. 

Observational evidence of this is not clear. Accord- 
ing to some studies, both in our galaxy and in others, 
the metallicity gradient measured using planetary neb- 



ulae (PN) is indeed the same (or even steeper) than 
the gradient measured with open clusters or HII regions 
jMaciel ot al' '2005'; 'StauEhcUini ot al.' I2OI0I : iMagrini et all 
|2009; Hornandcz-Martmoz ot al. 2009). According to other 
studies, however, the PN metallicity gradient is shallowe r 
than that of HII regions (|Perinotto fc MorbidellH bOOd V 
IYu et all (I2OI2I ) found that the gradient steepens with the ef- 
fective temperature of stars that have been used to measure 
the gradient. This again indicates that the gradient becomes 
shallower when we consider older objects, because of radial 
migration. 

Thus, to have the PAH content correlated with the local 
ISM metallicity in the "AGB" scenario, we need three con- 
ditions to be satisfied. First, the efficiency of the PAH syn- 
thesis should be related to the metallicity of the parent AGB 
star. Second, the parent AGB star should stay in approxi- 
mately the same metallicity environment during its lifetime. 
Third, PAHs should not leave this environment due to radial 
gas flows. Note that azimuthal mixing may also be impor- 
tant in this respect, as at least in the Milky Way there a re 
some azimuthal metallicity variations (|Balser et al.ll201ll ). 

There is one more argument which may turn out to 
be contradictory to the "AGB" scenario. Supposing that 
VSGs are represented by PAH clusters, it is reasonable to 
assume that their formation processes are similar to those 
of PAHs. Thus, we may expect a correlation between P24 
and metallicity, similar to what we see between Pg.o and 
metallicity. However, there are no signs for such correlation 
(Fig.Et). 

This argument can be irrelevant if at some circum- 
stances large grains over-power very small grains at 24 
\im. However, most studies where distinction is made be- 
tween grains of various types indicate that in the gen- 
eral ISM most significant contribution at 24 \im comes 
from VSGs (iDesert et al.lll990l : ICompie gnc ct al 2011). The 
same seems to be true in PDRs as well (Bcrnc ct al. 20071 ; 
ICompiegne et aLll2008l ) . In the immedia te vicinity of HII re - 
gions the situation is more complicated. iFlagev et all (|201ll ) 
have shown that under certain assumption on the grain size 
distributions in some regions within the Eagle Nebula emis- 
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sion from big grains (BG) i s indeed comparable o r even 
exceeds emission from VSGs. iDeharveng et alj l|2010l ) anal- 
ysed emission maps of several HII regions at 8 and 24 \im 
and stated that the source of emission at 24 mum (BGs vs 
VSGs) cannot be indentified clearly. We have initiated our 
own modelling of the infrared emission from an HII region 
that takes into account grain destruction and found that the 
contribution from big grains at 24 yuai is small in comparison 
to the VSG contribution (M. Kirsanova et al., submitted). 
Still, the question of the 24 [im emission source is far from 
being solved. 

In the "destructive" scenario, there is only one condi- 
tion that needs to be met. Namely, the PAH destruction 
efficiency should be related to the local metallicity, i.e., 
through the intensity and/or hardness of the UV field. The 
relation between qpah and the UV intensity (characterized 
by (7min) does seem to exist, albeit with significant scatter 
(Fig. [5^). Regarding on the scatter in the qpah — (O/H) and 
qpAH — Uuiin relation, we must note that there are two un- 
knowns in the "destructive" scenario that are able to smear 
out metallicity trends. These are the age of the HII region, 
and the initial PAH content qp^n, prior to the onset of star 
formation in the region. 

The absence of the P24 correlation with (O/H) (Fig. O 
does not confiict with the "destructive" scenario (if the emis- 
sion at 24 (j.m is indeed dominated by hot, small grains). 
It can be explained by the greater stability of small grains 
against the ultraviolet radiation field, compared to PAHs. 

Some arguments in favour of the "destructive" scenario 
are presented in the spatially resolved mid-IR study of the 
Small Magellanic Cloud (SMC) by Bolatto ct al. (2Q03) and 
ISandstrom et al.l (|2010l ). These authors have found a large 
range of the F8.0/-F24 values across the SMC star-forming 
regions despite their relatively uniform metallicity and con- 
cluded that the metallicity-dependent formation cannot be 
the pr imary factor, det e rmini ng the PAH abundance in the 
SMC. ISandstrom et all (|2010l ') also found that the FS.0/F24 
ratio correlates with gpAH, albeit with significant scatter. 
In our study the correlation appears to be stronger, but we 
only consider star-forming complexes, while in the SMC the 
greatest scatter in FS.0/F24 is ap parently observed in th e 
diffuse medium with gpAH < 1.5% (ISandstrom et al.ll2010l ). 

The factors which determine the initial PAH con- 
tent in an HII region are very diverse. First, the initial 
PAH molecules in a molecular cloud may be provided by 
AGB stars. In this case, the possible pre-existing gpAH- 
(O/H) correlation, established globally due to metallicity- 
dependent PAH formation, can be further enhanced locally 
due to metallicity-dependent PAH destruction. Second, even 
if PAHs synthesized in AGB stars are completely destroyed 
as they travel through the intercloud space ((?pah — 0), 
they may be synthesized in molecular clouds. As far as 
we know, there are no studies that address the forma- 
tion of PAHs in molecular clo uds directly. However, th ere 
are other related studies, like iBettens fc HerbstI l| 19971 ). If 
PAHs are indeed able to form in molecular clouds, their ob- 
served content is determined by the balance between the lo- 
cal formation (possibly metallicity-dependent) and the local 
(metallicity-dependent) destruction. Third, PAHs observed 
in star-forming regions may represent an intermediate stage 
in the metallicity-dependent destruction of larger particles 



(e.g., graphitic very small grains or amorphous hydrocarbon 
conglomerations) . 

To clarify the relative importance of all these factors, 
a detailed numerical model is needed (D. Wiebe et al., in 
preparation). Also, it would be useful to find a way to esti- 
mate the VSG mass fraction in HII complexes, and to ex- 
plore whether or not it shows similar correlations as the PAH 
mass fraction. 



6 SUMMARY AND CONCLUSIONS 

In this work we performed aperture photometry for more 
than 200 HII complexes in 24 spatially-resolved nearby 
galaxies from the SINGS sample. Spitzer and Herschel data 
are used to estimate the PAH abundances and the UV field 
intensity in these HII complexes. Using literature values of 
the metallicities of the HII complexes, we show that the well- 
known correlation between the PAH content and metallic- 
ity previously found for entire galaxies is preserved also for 
smaller-structure components, like star-formation complexes 
and giant HII regions. We lack galaxies with 12 + (O/H) 
around the threshold metallicity in our sample, so we are 
not able to trace the region in the qpah — (O/H) parameter 
space where qpah begins to increase with metallicity. Several 
arguments are presented in favour of the "destructive" sce- 
nario of PAH evolution, which relates the lower PAH mass 
fraction in low-metallicity environments to a stronger UV 
field. However, many factors are involved in PAHs evolu- 
tion, so more detailed modelling should be able to provide 
more clear information on their evolutionary path. 
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Table 2: Parameters of HII complexes. Pull table is available online. 



No. Object 


q(J2000) 


(5(J2000) 


aperture 


F3.6, mJy 


F4.5, mJy 


F5.8, mJy 


Fg.o, mJy 


F24, mJy 


F70, mJy 


F160, mJy 






h m s 




radius, " 
















IC2574 


1 


HI 


10 28 55.60 


+68 27 54.8 


12.0 


0.01 ± 0.01 


0.01 ± 0.01 


0.07 ± 0.01 


0.15 ± 0.01 


2.3 ± 0.1 


3.6 ± 1.8 


2.2 ± 2.0 


2 


H2 


10 28 58.91 


+68 28 27.7 


14.0 


0.12 ± 0.01 


0.10 ± 0.01 


0.21 ± 0.02 


0.31 ± 0.01 


2.7 ± 0.1 


45.9 ± 2.4 


55.1 ± 6.1 


3 


H3 


10 28 48.50 


+68 28 02.3 


13.0 


0.34 ± 0.03 


0.36 ± 0.02 


1.10 ± 0.03 


2.4 ± 0.1 


27.2 ± 0.3 


259.4 ± 6.3 


225.9 ± 8.5 


4 


H5-6 


10 28 50.18 


+68 28 23.3 


14.0 


0.19 ± 0.03 


0.17 ± 0.02 


0.33 ± 0.04 


0.63 ± 0.07 


5.3 ± 0.6 


129.7 ± 6.6 


111.3 ± 6.8 


5 


H8 


10 28 43.65 


+68 28 26.4 


12.0 


0.22 ± 0.04 


0.29 ± 0.04 


0.30 ± 0.04 


0.51 ± 0.03 


8.0 ± 0.2 


60.7 ± 4.0 


35.1 ± 4.7 


6 


HIO 


10 28 39.30 


+68 28 06.9 


12.0 


0.16 ± 0.02 


0.16 ± 0.02 


0.19 ± 0.01 


0.31 ± 0.01 


2.7 ± 0.2 


44.5 ± 3.6 


43.2 ± 3.4 


7 


H13-14 


10 28 30.84 


+68 28 08.4 


12.0 


0.06 ± 0.01 


0.04 ± 0.01 


0.08 ± 0.01 


0.09 ± 0.01 


0.8 ± 0.1 


8.7 ± 0.8 


11.3 ± 2.8 


8 


II 


10 28 48.40 


+68 28 02.0 


14.0 


0.35 ± 0.03 


0.36 ± 0.02 


1.10 ± 0.03 


2.5 ± 0.1 


28.0 ± 0.3 


273.1 ± 7.2 


236.7 ± 9.2 


9 


III 


10 28 50.91 


+68 25 26.1 


14.0 


0.34 ± 0.02 


0.30 ± 0.01 


0.52 ± 0.24 


0.84 ± 0.03 


4.1 ± 0.2 


104.2 ± 2.9 


112.1 ± 4.2 
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Table 3: Fitting parameters of HII complexes and rnctallicitcs. Full table is available online. 

No.Objoct gr-Aii 1% f/min 12 + log(0/H) 12 + log(0/H) 

KKO i PTO.") 



IC2574 



1 


HI 




30.0+i!'; 




8.14 


± 


0.04 


7.71 ± 0.05 


2 


H2 


0.5irr 


1 o + O.l 


le.o+r* 


8.09 


± 


0.04 


7.72 ± 0.05 


3 


H3 


"J- ' -0.1 


^•■-'-0.1 


9.6±oi 


8.15 


± 


0.04 


7.72 ± 0.05 


4 


H5-6 


0.5irr 




4.8t;;:i 


8.21 


± 


0.04 


7.82 ± 0.06 


5 


H8 




8.4i;];| 


i2.oii;i 


8.14 


± 


0.04 


7.72 ± 0.05 


6 


HIO 




1 q + O.O 




8.11 


± 


0.04 


7.70 ± 0.05 


7 


H13-14 




^•O-O.l 


^•^-0.2 


8.09 


± 


0.03 


7.66 ± 0.05 


8 


II 




^•O-O.l 


8.0«-J 


8.27 


± 


0.02 


7.97 ± 0.03 


9 


III 


o.7l°:l 


1 3+°-i 




8.35 


± 


0.03 


8.11 ±0.05 



